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^^ . The first lattice QCD numerical study of heavy quark-quark potentials at finite temper- 

t*^ ' ature is reported. Using the quenched approximation, we evaluate the color anti-symmetric 

("V^ ' and symmetric potentials. 

(N 

^ ■ Recently, the diquark (a two-quark system) has attracted much attention with 

regard to high-energy phenomenology. Jaffe and Wilczek proposed that the recently 
\^ \ discovered penta-quark state {0~^Y' is a bound state of {ud){ud)s, where (ud) rep- 

^^ ' resents for highly correlated u and d quark pairs. ^^ It is believed that at high baryon 

number density and low temperature there exists a family of color superconducting 

^ I phases, due to the quark pairing driven by the BCS mechanism.^) ( See Ref. 4) for 

^^ ■ a review of the history of diquarks and their role in high-energy reactions. ) 

(<— ^ . A quark-quark system is a color anti-triplet (anti-symmetric) or sextet (sym- 

VO ' metric) state: 

O ' 
Tt' 3x3 = 3* + 6, (1 



nxn = p + iin . 



B- 



o 
-)— » 

We believe that the quark-quark interaction is attractive and strong in the color 
jil^l anti-triplet channel, based on results obtained from the perturbation^^ and instan- 

^ ■ ton induced models. ^^ It is important to investigate the quark-quark potential using 

lattice QCD, which provides a non-perturbative and first principle basis for explor- 
ing the quark-quark interaction. To our knowledge, there has been only one such 
^ ' study, by Wetzorke and collaborators,"^' in which diquark correlation functions were 

H . calculated. 

In this work, we study the heavy quark-quark free energies at finite temper- 
ature that are obtained from the Polyakov line correlation (PLC). PLC was first 
investigated by McLerran and Svetitsky.^^ The free energy F is given by 



e-^^ 



Y.{<l>\e~^''\<l>), (2) 



where {(p) represents a state of gluons and heavy quarks. For the heavy quark-quark 
system of the color anti-triplet, we have 

\4>) = €abc^\xi,t = Q)^r{S2,t = 0)t|Gluons). (3) 

Here a, h and c are the color indices. Now, the summation in Eq. ((2)) is taken over 
a and all gluonic states. Using the relations^) 

V'(x,t) = LV'(^,0), (4) 
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{r{x, f = 0), V'(y, t = 0)t} = 6a,bd^,y, (5) 



where 

Nt 



L{x) = llUo{x,t) (6) 

and 



t=i 



we obtain the free energies for the qq sector in the symmetric and anti-symmetric 
channels, 

exp{-PFsiR)) = ^{TtL{xi)TiL{x2)) + ^{TiL{xi)L{x2)), 
exp{-PFas{R)) = -(TVL(xi)T>L(x2))--(TVL(xi)L(x2)), 

where Fg and Fas are the differences in the free energy for the cases with and without 
the static color symmetric (anti-symmetric) quarks located at xi and X2, R = |xi— X2I 
and f3 = 1/T. 

The formula (jSJ is gauge dependent, and therefore it requires the gauge fixing. 
We employ the stochastic gauge fixing quantization (SGFQ) with the Lorentz-type 
gauge fixing proposed by Zwanziger:^' 

Here, a, D°^ and r/" correspond to a gauge parameter, a covariant derivative and 
Gaussian random noise, respectively. The lattice formulation and a more detailed 
explanation of this algorithm can be found in Refs. 10)-12). 

In the quenched approximation, the system possesses Z3 symmetry. The gauge 
action is invariant under the transformation U^lx) — > zU4{x) or z'^Ui^x) on a time- 
slice hyperplane, where z = exp(27ri/3). In this transformation, (TrL(xi)TrL(x2)) 
and {TtL{xi)L{x2)) obtain a factor of z'^ or z'^. After an infinite number of sweeps, 
the quantities in Eq. Q vanish because 1 + z^ + z^ = 0. This is true also for the 
Polyakov line expectation value itself, and often |TrL| or (Re(TrL)'^)^''^ is taken. 
The latter procedure brings the Polyakov lines around ib27r/3 into regions around 
the positive real axis. This may be justified because this area is realized in full 
QCD, due to the dynamical quark effect. Therefore we restrict our simulations to 
the region satisfying — 7r/3 < Phase L < vr/S. 

In order to calculate the color dependent Polyakov correlation functions at finite 
temperature, we performed the quenched SU{3) lattice simulation with the standard 
plaquette action. The values of the simulation parameters used are the same as in 
the previous study for the measurement of qq PLC functions in the QGP phase. ^'^^ 
The spatial lattice volume was 24'^, and the temporal lattice size was set to Nt = 6, 
which determines the temperature as T = 1/Nta, where a is the lattice spacing. 
The corresponding critical temperature is estimated to be Tc ~ 256 MeV in Ref. 
14). The system temperature is changed when we vary the lattice cutoff, which 
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Fig. 1. The behavior of qq free energies in the symmetric and anti-symmetric channels with 
Langevin step width At = 0.03. 



is determined by Monte Carlo renormalization group analysis. ^^^ All of the PLCs 
are measured every ten Langevin steps and normalized by the value (TrL(O))^. We 
obtained between 3000 to 10000 data points after approximately 3000 steps were 
discarded as thermalization. 

Typical behavior of the symmetric and anti-symmetric free energies at T/Tc = 
2.02, 3.04, 5.61 is displayed in Fig. ^ The symmetric channel gives a repulsive force, 
while the anti-symmetric one gives an attractive force. As the system temperature 
is varied, all potentials change but their variations are observed to be small. 

If there is a strong attractive force between quarks even at high temperature, 
diquark degrees of freedom should be taken into account for analyzing hadron produc- 
tion in high energy heavy ion collisions at RHIC and LHC. Therefore it is important 
to investigate the diquark force by the lattice QCD simulation. As a first step for 
such a study, we define the effective force 



fiR,T) 
2^2 



^{F{R + l)-F{R)). 



(10) 



The temperature dependence of the strength of this effective force is shown in Fig. 
121 As the temperature increases, this strength decreases in both channels. 

In the leading order perturbation (LOP), the coefficients of the color exchange 
terms in the symmetric and anti-symmetric channels are Cqq[6] = +5 and Cgq[3*] = 
— 3 respectively. We compare the force strengths in the symmetric and anti-symmetric 
channels and define their ratio Fs/Fas, which approaches the value Cqq[Q\/Cqq[2>*] = 
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Fig. 2. The temperature dependence of the effective force in the symmetric and anti-symmetric 
channels with Langevin step width At — 0.03. 



—0.5 at short distances. Figure displays values of this ratio evaluated at T/Tc = 
2.02, 3.04 and 5.61. We find that as the temperature increases, three values at short 
distances become comparable with the expected values: However, as the temper- 
ature decreases, it deviates from —0.5. There are several possible sources for this 
deviation at shorter distances for small T. One is numerical instability due to a long 
autocorrelation near Tc. If this is indeed the cause, then larger scale simulations in 
the vicinity of Tc are required. Another possibility is the renormalization of PLC, 
i.e., {TrL{R)L{0)) = Z exp{—F{R)/T). If Z differs significantly from 1, we should 
take it into account when comparing F with that obtained from the perturbation. 
Recently, the Bielefeld group has proposed that the proper normalization of PLC is 
determined by the comparison of the T = qq potential data with the corresponding 
T ^ qq singlet free energy at short distances. ^^^ However, the determination for 
the normalization of the qq channels using this method is not performed. 

In the SGFQ algorithm with a Lorentz-type gauge fixing term, the qq symmetric 
and anti-symmetric channels are studied by calculating the PLC functions. However, 
they are not gauge invariant. Figure 0] displays the gauge parameter a dependence 
of Fg and Fas ^or three values of a from 0.6 to 1.3. We find that the dependence on 
a in this range is small, and in particular, the basic features of the curves do not 
change. 

In this paper we have reported the first non-perturbative lattice QCD simulation 
of the qq symmetric and anti-symmetric free energies. Calculations were carried out 
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Fig. 3. The temperature dependence of the ratio Fs/ Fas with Langevin step width At — 0.03. The 
dash-dot line stands for the expected ratio, C'qq[6]/Cqq[3*] = —0.5 at short distances. 



for both gauge dependent channels in the framework of the stochastic quantization 
with a Lorentz-type gauge fixing term. We find that the anti-symmetric free en- 
ergy constitutes an attractive force, while the symmetric free energy constitutes an 
repulsive force. 
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Fig. 4. Gauge parameter a dependence of the free energies in the symmetric and anti-symmetric 
channels at T/n = 3.04 {At = 0.03 ). 
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